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We consider a typical setup of cavity QED consisting of a two-level atom 
interacting strongly with a single resonant electromagnetic field mode inside a 
cavity. The cavity is resonantly driven and the output undergoes continuous 
homodyne measurements. We derive an explicit expression for the state of 
the system conditional on a discrete photocount record. This expression takes 
a particularly simple form if the system is initially in the steady state. As 
. a byproduct, we derive a general formula for the steady state that had been 

CN , conjectured before in the strong driving limit. 

: 

s . 

> ■ I. INTRODUCTION 

m : 

Recently there has been much experimental progress in sing le-atom cavity QED mm. In 
addition to their inherent fundamental importance, these experiments provide insight into 
the physics of open quantum systems, with potential applications to, e.g., quantum chaos ||, 
; quantum control ||, and quantum computing [§]. 

In this article we consider a typical experimental setup of single-atom cavity QED ||, 
as illustrated in Fig. 111. The setup consists of a single two-level atom located inside a high 
finesse optical cavity, which is externally driven. A set of photodetectors is arranged to 
monitor the field escaping from the system into the environment. We assume that the 
leakage of photons from the cavity mode through an output mirror is the only significant 
channel through which the system interacts with the environment. This assumption can be 
very realistic for high-finesse cavities 0. Also, for simplicity, we adjust the cavity length and 
the frequency of the driving field so that they both coincide with the frequency of the atomic 
transition. The cavity output is monitored using continuous homodyne measurements ||. 
These measurements are parameterized by one complex parameter: the reference field (3 
^ ■ which is added to the cavity output on a beam-splitter prior to the detection. 

Given the output of the photodetectors, it is possible in principle to write the conditional 
quantum state inside the cavity as a function of time and the measurement record. Usually, 
the conditional state is computed numerically 0,|TIJ using the formalism of stochastic mas- 



ter equations |pT| , p!2| ; these numerical computations can require very large computational 
resources. For some experiments, however, the ability to process data in real time is cru- 
cial H]. It is therefore important to develop analytical tools for conditional state evolution. 

In this paper we derive explicit expressions for the state conditioned on a discrete homo- 
dyne measurement record in the strong coupling regime, where the atom is strongly coupled 
to the intracavity field. Our calculations are valid for arbitrary driving field strengths. Our 
results are applicable to experiments such as the atomic cavity microscope P],0], where the 
strong coupling is essential, but strong driving leads to the problem of saturation J7|. We 
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give special attention to the experimentally important case that the system is initially in 
the steady state. 

The paper is organized as follows. In Sec. [TI| we describe the equations that model the 
physical system. In Sec. |TTT| we review the formalism of conditional quantum evolution and 
introduce the approximations for the strong coupling regime. In Sec. [TV] we derive a general 
formula for the system state conditioned on a discrete photocount record, for an arbitrary 
initial system state. In Sec. [V], we give a derivation of a general expression for the steady 
state. In Sec. [VI] we find simple formulas for the conditional evolution in the case that the 
system is initially in the steady state. We conclude in Sec. |V1J] , 




FIG. 1. Homodyne measurements in cavity QED. Basic parameters of the system are the 
strength of the atom-cavity coupling g, and the cavity field decay rate 7. The cavity is reso- 
nantly driven by an external laser field E, and the cavity output field is analyzed by the detectors 
D\ and D2 after being added to the reference field (5 on the beam-splitter. 



II. MATHEMATICAL MODEL AND MAIN APPROXIMATIONS 

Let |g) and |e) be the ground and excited states of the atom. For simplicity, we choose 
the cavity length so that the frequency of the resonant optical mode coincides with the 
frequency of the atomic transition. Using the dipole and the rotating wave approximations 
the interaction of the two-level atom with the electromagnetic field inside the cavity is 



described by the Hamiltonian [13 



H int = ig{a ] a - aa ] ) , (1) 

where a = |g)(e|, g is the strength of the atom-cavity coupling, and a is the annihilation 
operator for the intracavity field. Including dissipation and on-resonant driving of the cavity 
mode, the total unconditional master equation in a frame rotating at the driving laser 
frequency reads 

p = [—iH- mt + E(a) — a), p] + —(2apa^ — a)ap — pa)a) , (2) 

where p is the joint density operator for the atom and the intracavity field, E is the strength 
of the driving, and 7 is the rate of energy loss due to the leakage of photons from the cavity 
mode through an output mirror. 

From the experimental point of view the question of the steady state is very important. 
In fact, using contemporary techniques it is very difficult to prepare the system in question 
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in any other state. Using the Jaynes-Cummings model, Alsing and Carmichael || have 
shown numerically that in the strong driving limit E 3> g the system approaches a steady 
state of the form 

Pss = ^(|ar; +)(a; +| + |a*; -)<a*; -|) , (3) 
where |a; +) and |a*; — ) are two orthogonal quantum states 

\a;+) = _| a )(|g) +i| e )) = |a)|+) , 

|a*;-> = -^|a*>(|g)-*|e» = |a>|-> l (4) 

and where |a) is the coherent field state with amplitude 

a = (2E + ig)/ 1 . (5) 

This result has been confirmed in a more recent numerical simulation flO|| . Using matrix 
notation for the intra-atomic degrees of freedom in the basis {|±)}, Eq. @ can be rewritten 
in the convenient form 



_ x \a/\a\ u , v 


















which will be useful below. 

In this paper, we work in the strong-coupling regime (g 3> 7), which justifies considering 
the evolution on time scales large compared to l/g. In Sec. we give an analytical proof 
that on those timescales, Eq. (H) is a steady state of Eq. (0), for arbitrary values of the 
driving i£. 



III. THE MEASUREMENT 

We now rewrite Eq. (^) in the form 

p = Cp, (7) 

where the superoperator C is defined as 

£p = [— ii?i nt + E(a) — a), p] + —(2apa) — a^ap — pa) a) . (8) 

Let the initial condition be p(0) = po- Given superoperators So, J\ and J7" 2 such that 

S Q (t) = eV-*- J &, (9) 

the solution to Eq. (0) can be written using a Dyson expansion, 

00 

p{At)=J2 E p(h,...,k m ;At)p c (k 1 ,...,k m ;At) , (10) 

m=0 ki,...,k m 
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where trp c (fci, . . . , k m ; At) = 1 and 



p(ki, ...,k m ; At)p c (h, ...,k m ; At) 

f*2 



/•At 




/ dt m 


■ r dt 2 f 


lo 


Jo Jo 



Following |pj^2| we define the "smooth evolution" operator So as 

S (t)p = N {t)p[N {t)]t , 

where 



N (t) = exp 



-%H mt t + E(a ] - a)t - J(a f a + \(3\ 2 )t 



and the "jump" operators J\ and J 2 as 

J k p = C kP Cl, where C k = e^ k -^ 2 [a + {-l) k (3} . 

The following lemma, included for completeness, shows that the definitions of Sq. 
J 2 just given are consistent with Eq. (Rp. 

Lemma 1 The above definitions satisfy the requirement 



(11) 
(12) 
(13) 

(14) 
Ji and 



S (t) 



(15) 



and therefore Eqns. and fill] ) indeed give a solution to (^j. 
Proof 

Keeping terms to first order in r we have 

S (t)p = N (t) P [N (t)^ 

= p + {[-iH int + E(a ] - a), p] 

- Itfap + pa^a) - i\P\ 2 P )t + 0(r 2 ) 

= (1 + r • C)p - l(apa^ + \/3\ 2 p)r + 0(r 2 ) . 

On the other hand by direct calculation we have 

J kP = l[apa} + (-l) fc (/W + (3*ap) + 2 p] 

which implies that 

(Ji + J 2 )p = l{apa) + \{3\ 2 p) . 
Equation (|16"D therefore becomes 

S (r)p = (1 + r • [C - (Ji + J 2 )])p + 0(r 2 ) 
Taking the limit r^Owe have Eq. fll5|) as required. □ 



(16) 



(17) 



(18) 



(19) 
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There are many different definitions of Sq, J\ and J2 that satisfy the above lemma. How- 
ever, definitions fll2|) and (14) are somewhat special: the quantities p c {k\, . . . , k m ; At) and 
p(ki, . . . , k m ; At) which they define have an important physical meaning [O. Suppose that 
the continuous measurements were performed over the time interval At and recorded as a se- 
quence (k\, . . . , k m ; At) of photodetector labels in the order of photodetections. For example, 
kj = 1 would mean that the jth photodetection was registered by the first detector. Then 
the probability of the measurement record (ki, . . . , k m ; At) is given by p(ki, . . . , k m ; At), and 
the corresponding conditional state is p c {k\, . . . , k m \ At). 

We will now prepare to consider the conditional system evolution on time scales large 
compared to 1/g. First we notice that 



where 



We define 



where 



H mt = H + H U (20) 



H = -g(a* + a)a y /2, a y = i{a^-a) > 2 _ 
Hi = ig(aJ — a)a x /2 , a x = + a . 



Q = exp(-W t - Wit + Ft) , (22) 



F = E(a^ - a) -ja^a/2. (23) 



These definitions are connected to the definition fll2"|) of the smooth evolution operator via 
the relation 



We rewrite Q in the form 



N = e~ 7l/3|2 * /2 Q . (24) 



Q = e - iHat R , (25) 



so that 



^ = -iHoe^Ro + e~ iH ^ . (26) 
dt dt v ; 



From the definition (g2|) and Eq. ( 25Q we have: 

= {F _ iHo _ W^e'^Ro . (27) 
Combining the last two equations we obtain that Rq obeys the equation 

^ = (X(t) + e iHot Fe- iHot ) R , (28) 

where 
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X(t) = -e tHot iH ie - iHot . (29) 

Using the Corollary to Theorem 1 from the Appendix together with the identity e~ ya 1 a = 
(a + y)e~ yat we obtain 

2X{t) = ge iHot {a) - a)a x e- iHot 

= g(a*-a- iga y t)e iHot a x e- iHot . (30) 

The identity e lAeTy = cos A + ia y sin A gives 

e iHot a x e- iHot = a x cos[gt(a^ + a)} - a z sin[^t(a t + a)} , (31) 

where a z = ia y a x . Finally we obtain 

2X(t) = g(a l — a — ig<r y t) (a x cos[gt(a^ + a)} — a z sm[gt(af + a)]) . (32) 



At time scales large compared to 1/g we can neglect terms oscillating at frequency 1/g in 
Eq. (p8[). This means we can set X(t) = in Eq. (p8|), which becomes 

^^e iHot Fe- iHot R . (33) 

This approximation has some similarity with the standard rotating-wave approximation. 
Now consider the operator 

M = exp(-iH t + Ft) . (34) 
Using the same technique as in Eqns. (|25| - 128|) it is easy to show that 

M = e~ iH ^R x , (35) 

where R\ obeys the equation: 

dRi 



dt 



e iH t Fe -iH Q t Ri _ (gg) 



This equation coincides with (^), which means that at time scales 5t ^> 1/g the operator 
Q can be replaced with M. The smooth evolution Sq can therefore be approximated as 

S ^S, (37) 

where S is defined as 

S(t)p = N(t)p[N{t)]i, (38) 

and where 



N(t) = exp 



-iH t + E(a ] - a)t - |(a f a + \f3\ 2 )t 



(39) 
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IV. CONDITIONAL EVOLUTION FOR ARBITRARY INITIAL STATES 



In this section we derive a general formula for the state, p c (k\, . . . , k m ; At), conditioned 
on a discrete photocount record for an arbitrary initial state. The formula is a direct 
consequence of two technical theorems, whose proofs are given in the Appendix. At time 
scales 5t ^> 1/g, the theorems allow us to simplify Eq. ([□]) by changing the order in which 
the smooth evolution operators S and the jump operators appear. 

Using these theorems we can proceed with the calculation of the conditional density 
matrix p c {k x , k m ; At). We have from Eqns. ([□]), ([TJ) and Eqns. (|37|-[3|) that 



p(ki, ...,k m ; At)p c (ki, ...,k m ; At) 

r At 



/■At 


rAt , 


1 dt m • 


■ dt 2 




Jo Jo 



iv(*2-fi)cktf(fi)W 



(40) 



We can now use Theorem 2 to compute the operator in the square brackets. We have, for 
instance, 



C fcl iV(*i) = N(t!)f k 



3 -7*l/2 



-(2E + iga y ) + (-l) k ^ 
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(41) 



Then, using the identity N(t 2 — ti)N(ti] 
of calculations we have 



-/V(i 2 ), we see that repeating the same type 



N{At - t m )C k 



N(At) n h v 

P =i 



N(t 2 -t 1 )C kl N(t 1 ) = 

I _ e -7tp/2 • 
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Using the identity fkf k = 7/2 we therefore have 



(42) 



p(ki 



k m ] At) p c {k\, . . . , k m \ At) 



r 



2 m m\ 



N(At)G(p ,P)N^(At) 



where 



G(p ,P) 



At 



dt r 



At 



n 

P =i 



-jt p /2 



dt m —i 
1 

a H 



At 

dti 

(l -7*p/2 
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(2E + tga y ) + (-l) k r(3 



Po 



(43) 



(44) 



For notational convenience, we do not indicate explicitely the dependence of G(p , (3) on the 
measurement record {k\, . . . , k m ; At) which, however, should always be remembered. 

Equations ( |43"|) and (0) have a relatively simple structure. The terms N(At), which are 
given in factored form by Theorem 1, are the same for all possible measurement records. This 
means that all the information about the measurement records is contained in the function 
G(po,j3). The integrand in G(po,f3) is a polynomial in a, a y and po- The scalar coefficients 
of this polynomial are constants or proportional to either e _7 * p / 2 or e~ ltp . Therefore all the 
integrals in Eq. ( f4"4| ) can be easily evaluated, so that G(p , (3) takes the form of a polynomial 
in a, o y and p with known coefficients. In this way, Eqns. ([43]) and (0) provide an explicit 
solution for the conditional evolution on the time scales considered. 
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V. DERIVATION OF THE STEADY STATE 



In this section we show that, at timescales 5t ^> l/g, the state p ss defined by (KJ) is a 
steady state of the master equation (|7]). Notice that the only free parameter in our homodyne 
measurements is the complex parameter (3. If we can find a value of /3 such that for any 
measurement record (k\, . . . , k m ; At), the conditional density matrix satisfies 

p c (kx, . . . , k m ; At) = p ss , (45) 

then p ss must be a steady state. This is because the solution ( ]I0D of the unconditional 
master equation ((7|) becomes, in this case, 

oo 

p(At)=J2 E p(fci,...,fc m ;At)p aa = p ss (46) 

m=0 fci,...,fc m 

for any At. Intuitively, one would expect that, if subjected to a nontrivial measurement, 
the system would normally depart from the steady state. In our case, however, we will find 
that Eq. (f45|) is satisfied for all real values of (3. 

Before we proceed with our rigorous analysis it may be helpful to develop some intuition 
about the dependence of the conditional evolution on (3. In particular we are interested in 
the dependence of the conditional evolution on the phase <fi = arg (3. In our analysis we deal 
with the conditional evolution conditioned on a discrete photocount record, which is the most 
general case. However, a lot of insight about the dependence of the conditional evolution 
on the phase <fi can be gained by taking the limit \(3\ — > oo. In this limit the detectors 
are registering continuous photocurrents rather than discrete photocounts. Because the 
resulting measurement records can be viewed as continuous functions of time it becomes 
possible to derive a master equation for the conditional density matrix p c . According to 
Ref. |TI|] , this can be done by taking the double limit \(3\ oc e _1 — > oo and 7 At oc e 3 / 2 — > in 



the Dyson expansion fllPl). If the measurement record consists of the difference photocurrent 
J_ = 1% — It, where I\ and I2 are the photocurrents detected by the first and the second 
detectors respectively, then the resulting master equation for the conditional density matrix 
becomes WWU 

p c = Cp + v^(e-^apc + e^p c « f - tr [p c (e"^a + e*V)] p c )£ , (47) 



where rj is the efficiency of the photodetection, and £ is the Gaussian white noise which, in 
practice, should be taken from experimental observations of the difference photocurrent I 
via the relation 



I_ = \{3\ (j V tr[p c (e^a t + e^a)] + . (48) 

Compared to the unconditional master equation (0), Eq. (O) has an additional term 



Vj[e-^ap c + e l<t> p c a) - tr p c (e"^a + e l V) p c )£ , (49) 



which, for = and p c = p ss , is proportional to p ss . This means that, if £p ss = 0, i.e., 
if pss is a steady state of the unconditional evolution, then conditional and unconditional 



S 



evolution coincide for = 0. This situation is similar to the one described by Eqns. fl4"5|) 
and (|46|) , which suggests to consider the case of real f3 in the following rigorous derivation. 

We now substitute p = p ss from Eq. (^) into Eq. (Q), keeping (3 arbitrary for the 
moment. We obtain 



G( Paa ,0) = dt m dtm-x-'-l dtJ]le-^ 2 [a + f(a y ,P)))pJ---y, (50) 

<^ <^ p=X 

where 

f(a y , (3) = e7W '~ 1 (2E + iga y ) + (-l) k ^e^ 2 . (51) 

We note that 

[a y , Pss } = and {a y f = 1 . (52) 

Using the first of these properties and the expression for p ss as given by Eq. @ , we have by 
direct calculation 



[a + f(a y , (3)}p ss [a + f(a y , /3)]t = [f(a yi (3) [f(a y , /3)]t + 2Re[f(a y , (3)}Re(a) 

+2lm[f(a y , /3)}lm(a)a y + \a\ 2 )p ss , (53) 



where a = (2E + ig)/^. We will use this equation for imaginary (3 in the next section. 
For the rest of this section, we assume that (3 is real. Using this and the fact that 



{(Ty) 2 = 1, we find that 



fMlfMV = AE2 ^ ^ [e 7W2 - l f + P 2 ^ tp + -^(-l) fcp (e^ - e^ /2 ) , (54) 



and 



Rei/KMRe,*) + Lo^W^, = **±A<*» - !) + ^(-1)^ . (55 ) 
Because | ct | 2 = \2E + ig\ 2 j^ 2 = (4E 2 + g 2 )/^ 2 we therefore have according to Eq. (|53|): 

[a + f(<r v , P)]p BB [a + f(a y , 0ft = e* ^* 9 " + (-l) fc ^ + /3 2 ] . (56) 
Substituting this into (ISO) we obtain: 



7 2 7 



G(p ss ,/?) = (At)- ft (^!±^! + (-i)^l^ + ^) pss . (57) 
P =i v 7 7 7 

Therefore, according to Eq. (|43[) , 

p c (k u ...,k m ;At) ocN(At)p ss N\At), for any real /3 . (58) 
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As the final step of our argument, we now prove a lemma that, together with Eq. (|58|) and 



the normalization of the density matrix, implies Eq. (45). 



Lemma 2 Smooth evolution leaves p ss invariant in the following sense: 

N{At)p SB N\At)ocp ss . (59) 

Proof 

Because a y and p ss commute, we can see from Eq. (^|) that the smooth evolution leaves p ss 
diagonal: 

N(At)p»[N(At)]i = ( A Q X I ) . (60) 
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Using Theorem 1 we have 



2e- 2Z ^A 1 = (e-^ a e z ^'e z ^)\a){a\{ 



2e- 2Z iA 2 =(e-^ a e z * a 'e z s a )\a*)(a*\(---)\ (61) 



where 



7 

Zf = H±ll( e -^t/2 _ x) ( g 2) 



In order to calculate A 1 we use the identity e Aat \a) (a|e A * a = e' Q+A ' 2 ' a ' 2 |a + A) (a + A| which 
gives 

2 e - 2 *A 1 = |e z 3-«| e l^ + l 2 H«l 2 (e-^ ata |a + Z 2 + )(a + Z 2 + |e- 2 ^ ata ) . (63) 
Now, with the help of the identity e~ Aata |a)(a|e _Aata = e' Q ' 2(e 2A ~ 1 - ) |ae~ A )(ae~ A | we have 

2e~ 2Z ^A 1 = |e z 3~-|el a+z 2 + l 2e - 7At H Q l 2 |(a + Z+)e- 2 # i )((« + Z+)e- 2 # £ |. (64) 
Using the definition of Z 2 and the value of a = (2E + ig)/ r y we see that 

(a + Z^e' 2 ^ = a. (65) 

Therefore 

2e- 2Zl A 1 = \e z * a \\a)(a\. (66) 
Repeating the same arguments for A 2 we have from Eq. (pl"|): 

2e~ 2Zl A 2 = \e z z a *\\a*)(a*\ . (67) 
Because \e Zs a \ = \e z s a *\ we can now see that 

(Ax \ / \a)(a\ \ , . 

( o aJ^V l«*)(«lj =2pss - (68) 

Together with Eq. (|60|) this completes the proof. □ 
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VI. CONDITIONAL EVOLUTION STARTING FROM THE STEADY STATE 



In the previous section we have shown that, for a real value of (3, a homodyne mea- 
surement does not give any information about the system once it has reached the steady 
state p ss . Although this fact was useful in confirming that p ss is indeed a steady state of the 
system, such a measurement would be pointless in practice. 

We therefore consider the case of purely imaginary (3, for which the homodyne measure- 
ment does provide information about the system. We write (3 in the form /3 = i(3$, where 
(3q is real. To find the conditional density matrix in this case, we go back to Eqns. (|50H53|) 
and obtain by direct calculation: 

fWVMf = ^^V W2 - l) 2 + fte* + (-l)*^(e* - e*>P) , (69) 



and 



Re[/K)]Re(a) + lm[f(a y )}lm(a)a v = ^ \ 9 \ e^ - 1) + (-l) k ^ e ^ 2 a y . (70) 



Therefore 

rAt t-At [At 

G(p ss ,i(3 )= dt m dt m _ x ---\ dti 
Jo Jo Jo 

ft (^±^ + Pi + (-l) fc ^[l + (a y - l)e-^/ 2 ]) Pss . (71) 
Performing the integration we obtain 

G( Pss , i/9b) = ft [At • {^±1- + (31 + (-1)^) 

p=i H H 

1_ e -7At/2 

+(-l) k Hg(3 (a y - 1) ] Pss . (72) 

Because p ss and a y are both diagonal in the basis {|±)} defined in Eq. (|J), the conditional 
density matrix can be written in the form 

^ A1 ) = ( A,l f wi)- < 73 > 

where a = (2E + ig)/j- For the eigenvalues Ai and A2 = 1 — Ai we have the following simple 
formula: 

*i = ff (74) 

where b = (4E 2 + g 2 + 7 2 /9q)/(2(77/9o)- Similarly simple expressions can be obtained for any 
complex reference field f3. 
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VII. CONCLUSION 



In this paper, we have given explicit formulas for the quantum state evolution conditioned 
on a discrete homodyne photocount record for a typical experimental setup in single-atom 
cavity QED. These formulas have potential applications for the real-time processing of ex- 
perimental data. The general methods developed here can be applied to a wide class of 
similar systems. For example, it should be straightforward to generalize our results to the 
case of heterodyne measurements. 
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APPENDIX 

Theorem 1 The operator 

M(t)=exp ig^-(a [ + a)t + E(a ] - a)t - ya f a (75) 
can be factorized as 

M(t) = e Zl e~^ a e Z2 ^e Z3a , (76) 

where 

^ = ^^(l-e-^- 7 t/2) 

z 2 = B±mL (e «/* _ 1) 

7. 

1) (77) 



2E-iga y /2 
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Proof 

Because a, a*, a^a and 1L span a Lie algebra, M(t) can be factorized in a systematic way as 
follows. First we find a function x(t) such that 

M(t) = e x{t)aU M{t) , (78) 

where M(t) is an exponential of a linear combination of a and . We will then repeat the 
same procedure factorizing M which will conclude the prove of the theorem. 
Equation (|7|) gives 

^^aW^M + e^f . (79) 
alt alt y ' 
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On the other hand, Eq. (ff5| ) gives 



< n[ [ig^a? + a)+E(a t - a) - ^a]e xafa M . (80) 



dt 2 v ' v '2 

Comparing this expression with the previous one we have 

dM 



dt 

where 



[ X (x)-(x + j/2)a^a]M , (81) 



x ( x ) = e ~^ a [ig^-{a) + a) + ^ - a)}e xa ' a . (82) 
Using the identity e - xa ^ a ae xa ^ a = ae x ) the above equation can be rewritten as 

X (x) = (£ + i^e" V — (E — ig^)e x a . (83) 
Looking at Eq. (]8T|) we demand that 

x + 7 /2 = , (84) 

thereby making dM/dt independent of a^a. From Eq. (|75|) we see that M(0) = 1 and 
therefore, we choose, in accordance with Eq. (|7S|), that 



x(0) = and M(0) = 1 . (85) 

With these conditions equation (Q) can be integrated to give, according to Eqns. ( [78] ) and 
(0), 

M(t) = e-^ ata M(t) , (86) 

where 

^ = [(£ + i<^)e^V _ ( E - %g °-JL)e'^ 2 a\e^ M' . (87) 

The proof of the theorem will be completed if we repeat the same procedure for factorizing 
M. As before we introduce a function y(t) such that 

M(t) = e y{t)at M'{t). (88) 

We therefore have 

d 4 = ya^M' + e>^. (89) 

dt y dt y 1 

Combined with Eq. fl87|) this gives 

—j— = [(E + ig^-)e <t/2 a) — (E — ig^e^e^ ae^ - ya^}M' . (90) 
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Using the identity e yat ae yat = a + y we rewrite the above expression as 



dM' 



[(E + ig^)^ 1 ' 2 - y]a^ — (E — ig^) e -^ 2 {a + y))M'. (91) 



dt V y "2' a> v 3 2 

We eliminate a) from this expression by setting 

y = (E + ig^' 2 . (92) 

Equation ( p8|) suggests the boundary conditions 

y(0) = and M'(0) = 1 . (93) 



Performing integration in (|9~1~D and in (|9~2"D using these boundary conditions and the fact that 
a 2 = 3L we have according to (|88|) 

M(t) = exp [^±^(e^ 2 - lWlM'm , (94) 

7 

where 

M'(t) = exp[ AE2 t 9 \ l - e-^ 2 - 7 */2)] exp[- 2E ~^ (l - e^/ 2 )a] . (95) 
7 2 '7 

This completes the proof of the theorem. □ 
Corollary 

e -iH t _ e -g 2 t 2 /8 e igta y a i ' /2 e igtcr y a/2 _ ^gg^j 

Proof 

This can be established easily by repeating the arguments of Theorem 1 for E = and 
7 = 0. □ 

Theorem 2 Using the definition 



h = ^e^ k -^ 2 (97) 
and the notation of Theorem 1, we have 

C k M(t) = M(t)f k [e-^ 2 a + (2E + iga y ) + (-l) k (3] . (98) 

7 

Proof 

By definition [Eqns. ( |i~4"D and (|97|) 1 and using Theorem 1 we have 



C*M (t) = f k e Zl [a + (-l) k p}e-^ a a e Z2a e Zsa , (99) 

where Zi, Z 2 and Z 3 are specified in the statement of Theorem 1. Using subsequently the 
identities e - xata ae xata = ae x anc [ then e~ yai ae yai = a + y we have 
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C k M(t) = f k e Zl e-^ ata [e^ t/2 a + (-l) k p]e z ^ e z * a 



f k e Zl e-^ a e Z2a \e^ t/2 (a + Z 2 ) + (-l) k f3}e Z3a 
= f k M(t)[e'^ 2 (a + Z 2 ) + {-\) k f3} . (100) 

Putting the value of Z 2 from Theorem 1 we have Eq. (P5|) as required. □ 
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